Kondo effect in a quantum critical ferromagnet 
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We study the Heisenberg ferromagnetic spin chain coupled with a boundary impurity. Via Bethe 
ansatz solution, it is found that (i) for J > 0, the impurity spin behaves as a diamagnetic center 
and is completely screened by 25* bulk spins in the ground state, no matter how large the impurity 
spin is; (ii) the specific heat of the local composite (impurity plus 25* bulk spins which form bound 
state with it) shows a simple power law C; oc ~ Ti; (iii)for J < 0, the impurity is locked into the 
critical behavior of the bulk. Possible phenomena in higher dimensions are discussed. 

c/3 i Kondo problem or the magnetic impurity problem in an electron host plays a very important role in modern 
. condensed matter physics. It represents a generic non-perturbationable example of the stppngly correlated many- 
body systems. Recently, with the development of research on some low-dimensional systemsu and the observation of 
unusual non-Fermi-liquid behavior in some Ueavy fermion compoundso, the interest in this problem has been largely 
■ renewed. The multi-channel Kondo aroblemS provided the first example of impurity systems whicb-.sliow non-Fermi- 
liquid behavior at low temperaturesu. In a Luttinger liquid, the impurity behaves rather diffcrenttrtl from that in a 
Fermi liquid and may interpolate between a local Fermi liquid ancLspme non-Fermi liquid.0 Some new quantum critical 
phenomena have also been predicted in some integrable modelstm. Generally speaking, these new findings indicate 
that the quantum impurity models renormalize to critical points corresponding to conformally invariant boundaxy 
conditional]. Another important progress is the study on the Kondo problem in Fermi systems with pseudo gapliil, 
i.e., the density of states p(e)-is power-law-dependent on the energy, p(e) ~ e r . With renormalization group (RG) 
["•"». ■ analysis, Withoff and FradkinHl showed that there is a critical value J c for the Kondo coupling constant J. For J > J c , 
\ Kondo effect occurs at low temperatures, while for J < J c , the impurity decouples from the host. We note that all 
• the quantum critical behaviors mentioned above only occur for T — > and therefore fall into the general category of 
quantum phase transitionsE3. ■— ■ 
In an earlier publication, Larkin and Mel'nikov studied the Kondo effect in an almost ferromagnetic metallic. 
With the traditional perturbation theory they showed that the impurity susceptibility is almost Curie type with 
logarithmic corrections at intermediately low temperatures. However, the critical behavior of a Kondo impurity in a 
quantum critical ferromagnet has never been touched. The main difficulty in approaching this problem is that almost 
all perturbation techniques fail in the critical regime and exact results are expected. As discussed in some recent 
worksEj, the critical behavior of the impurity strongly depends on the host properties and seems to be non-universal. 
Typical quantum critical ferromagnet is the Heisenberg system in reduced dimensions (d < 2). These systems have 
long-range-ordered ground states but are disorder at any finite temperatures due to the strong quantum fluctuations. 
In this paper, we study the critical behavior of an impurity spin coupled with a Heisenberg ferromagnetic chain. The 
model Hamiltonian we shall consider reads 



1 lf ~ 1 

H = -- V a 3 -.- 3j+i + Jai ■ S, (1) 



a 

where ffj is the Pauli matrices on site j; N is the length of the chain; S is the impurity spin sited at one end of the 
chain; J is a real constant which describes the Kondo coupling between the impurity and the host. The problem is 
interesting because (i)the model is not conformally invariant due to the nonlinear dispersion relation of the low-lying 
excitations, e(k) ~ k 2 , and p(e) ~ e _ 2 ; and represents a typical quantum critical system beyond the universality 
of the conventional Luttinger liquid; (ii)thje Hamiltonian is very simple (without any superfluous term) and allows 
exact solution via algebraic Bethe ansatztLx In fact, most known methodsEHa developed for the impurity problem in 
a Luttinger liquid can not be used for the present system due to the strong quantum fluctuations. 

Let us first summerize the solution of (1). Define the Lax operator Lj T (X) = X + i/2(l + Bj ■ r), where f is the 
Pauli matrices acting on the auxiliary space and A is the so-called spectral parameter. For the impurity, we define 
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Lo r = A + i(l/2 + 5 ■ t). Obviously, Lj T and Lq t satisfy the Yang-Baxter equation (YBE)IIa. It can be easily shown 
that the doubled-monodromy matrix 



T T (A) = L Nt (X) ■ ■ ■ L 1t (X)L 0t (X - ic)L 0T (X + ic)L lT (X) ■ ■ ■ L Nt (X), 
satisfies the reflection equation 

L TT >(X - fi)T T (X)L TT/ (X + fi)T T r(fj,) = T T ,(ji)L TT >(\ + u)T T (X)L TT ,(X - /i). 



(2) 



(3) 



From the above equation we can show that the transfer matrices (9(A) = Tr T T T (X) with different spectral parameters 
are commutative, [0(A), 0(fi)] = 0. Therefore, 0(X) serves as a generator of a variety of conserved quantities. The 
Hamiltonian Eq.(l) is given by 



h = l -J(-i) N -^e(X)\ x=0 + i(jv + i- J), 



with J = l/[c 2 — (5 + 1/2) 2 ]. Following the standard methoaQi we obtain the Bethe ansatz equation (BAE) 



2iV 



X 3 ' ~ h \ X 3 ~ + c ) X 3 ~ — c ) TT X i — Xi — i Xj + Xi 



X 3 + l) X 3 + { ( S + c ) X 3 +i(S -c) f-J: Xj -Xt+iXj + Xi+V 



n 



(4) 



(5) 



with the eigenvalue of Eq.(l) as 



M 



£({A,})=E 



\2 

A 3 



-(N-1) + JS, 



(G) 



where Xj represent the rapidities of the magnons and M the number of the magnons. _ 
Ground state. In the thermodynamic limit, the bulk solutions of Xj are described by the so-called n-stringsEZl. 
However, due to the presence of the impurity some boundary bound states may exist for c > S, which are usually 
called the n — A:-stringst3: 



XI 



S) 



k,k+ !,■ 



(7) 



In the ground state, only some n — 0-strings may survive. We call them boundary n-strings. In our case, n > has 
also an upper bound n < 2S — 1 since Xj = ±i(c + S) are forbidden as we can see from Eq.(5). No bulk strings 
can exist at zero temperature since they carry positive energy. Boundary bound state can exist only for c > S + 1/2 
(antiferromagnetic Kondo coupling) because in this case, the boundary n-strings carry negative energy. For zero 
external magnetic field, the most stable boundary string has the length of 2S with the energy 62s = 2S/[S 2 — (c— 1/2) 2 ]. 
Therefore, the impurity contributes a magnetization of — S. Such a phenomenon can be understood in a simple 
picture. Due to the antiferromagnetic coupling between the impurity and the bulk, 25 bulk spins are swallowed by 
the impurity at zero temperature to form a 25 + 1-body singlet. This singlet does not contribute to the magnetization 
of the ground state. In this sense, the impurity is completely screened, no matter how large the impurity moment is. 
Such a situation is very different from that of the conventional Kondo problem, where the impurity moment can only 
be partially screened by the host when 5 > 5' (5' the spin of the host particles)Eil This difference is certainly due 
to the different properties of the hosts. In the antiferromagnetic spin chain or a normal metal, the spin correlation of 
the bulk is antiferromagnetic type which repels more than one bulk spin or electron to screen the impurity. However, 
in a ferromagnetic spin chain, the bulk correlation is ferromagnetic which allows and in fact enhances some bulk spins 
to form a larger moment to screen the impurity. The local singlet is nothing but a bound state of 25 magnons. The 
boundary string may be broken by the external field. In fact, there are 25 critical fields 



m = 



25 



25- 



(c-i) 2 -5 2 (c-^±i) 5 



(5- §) 2 



n = 1,2, 



•,25. 



(8) 



When H r c l < H < H 7 c l+1 , only a boundary (25 — n)-string survives in the ground state and when H > i? 2S , any 
boundary string becomes unstable. Notice that at H = H™, the ground-state-magnetization has a jump SM = 1, which 
corresponds to some type of quantum phase transition. The finite value of Hi indicates that the zero temperature 
susceptibility of the local singlet is exactly zero. 
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Thermal BAE. Since we are interested mostly in the critical behavior, we consider T,H << U\ and J > case 
in the following text. In this case, any excitations breaking the boundary string can be plausibly omitted due to the 
energy gap associated with them. With the standard thermal Bethe ansatzO, we derive the thermal BAE as 



ln(l + ?7„) = 



27ra n (A) + nH 
T 



J2 A^Hl + ty-^A)], 



(9) 



rn—1 



or equivalently, 



lnr 7l (A) = - 5 (A) + Gln[l 



ln?7„(A) = G{ln[l + r?„+i(A)] + ln[l + »fc_i(A)]}, n > 1, 



(10) 



lim 



In rj n 



H „ 
= — = 2x , 



where a„(A) = n/27r[A 2 + (n/2) 2 ], A mn = [m + n] + 2[m + n - 2] H h 2[|m - n\ + 2] + [\m - n\}: g(X) = 1/2 cosh(7rA); 

r) n (\) are some functions which determine the free energy of the system; and [n] and G are integral operators with 
the kernels a n (A) and 5(A), respectively. The free energy is given by 



F bulk =F -(N + ~)T f <?(A){ln[l + Vl (X)} - + H }d\, 



(11) 



1 00 f 

Fim P = 2 T Ey K(A)ln[l + 7 7i ; 1 (A)]dA, 



where a„, m (A) = ^™ l "( m '™) a n+m+1 _ 2i (A); 4>' n (X) = a ni2 s(A - ic + i) + a„. 2 s(A + ic - i)\ F is the ground state 
energy; F^uik and Fi mp are the free energies of the bulk (including the bare boundary) and the impuritVp-jjespectively. 
Notice that Eq.(9) and Eq.(lO) are more difficult to handle than those o-£i^ ne antiferromagnetic chaintH, since here 
all tjn diverge as for T — ► 0. These equations were solved numerically! 2 !] in studying the critical behavior of the 
ferromagnetic Heisenherg chain. In addition, Schlottmann gave an analytical result based on a simple correlation- 
length approximation! 2 .!! and the result coincides with the numerical ones very well. As we can see from Eq.(9) and 
Eq.(10), when T — > 0, r) n — > 00. To arrive at the asymptotic solutions of r] n (X), we make the ansatz rj n {\) = 
exp[27ra n (A)/T](/) n . Substituting this ansatz into Eq.(10) we readily obtain (f> n ~ 1 for finite n and A. Therefore, 



i] n sa exp[ 



27ra»(A) 
T 1 



(12) 



On the other hand, when A — > 00 o: 
solution of rj n for very large A or n\ 



00, the driving term in Eq.(10) tends to zero. This gives another asymptotic 



Vn 



sinh 2 [(n + 1)ieo] 
sinh 2 xo 



- 1 



O(-e-^l), 



(13) 



For intermediate A and n we have a crossover regime. We call Eq.(12) as the strong-coupling solution, while Eq.(13) 
as the weak-coupling solution. By equating them we obtain two types of crossover scales, A c (n) for small n and n c (T), 



A c (n) 



4Tln(l 



n c {T) 



1 



4Tln(l 



AT In T ' 



(14) 



which characterize the crossover of the strong-coupling regime and the weak-coupling regime. Notice that the strong- 
coupling solution gives the corxact ground state energy and the low-temperature thermodynamics is mainly dominated 
by the weak-coupling solutiono. With such an approximation, the recursion for rj n can be performed by substituting 
the asymptotic solutions into the right hand side of Eq.(9) and therefore the leading order correction upon the 
asymptotic solutions can be obtained. In the following recursion process, we adopt the strong-coupling solution in 
the region of A < A c and n < ri c , while the weak-coupling solution is adopted in other cases. This corresponds to an 
abrupt crossover, which does not affect the temperature dependence of the thermodynamic quantities in leading orders 
but their amplitudes. For convenience, we define (n(A) = hi[l + 7?n(A)] — [27ra„(A) + nH]/T, which are responsible for 
the temperature-dependent part of the free energy. 
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Low-temperature susceptibility of the impurity. For convenience, we consider 2c = integer case. Taking the 
boundary string into account, the free energy of the impurity can be rewritten as 



F imp = \t J .9(A)[C 2c+ 25- 2 (A) - sgn(2c~ 25 - 2)C| 2c -2S-2|(A)]dA. (15) 

Substituting the asymptotic solutions Eq.(12) and Eq.(13) into Eq.(9) and omitting the exponentially small terms, 
we obtain 

n c -. rj />OC p— A c (m) 

C»(A)^^{ln[l+ ]--^}[/ + / ]A mn (X - X')dX' 

sinh(l + n c )x 

+2n c ln — , (16 

Sinn n c XQ 

where A mn is the kernel of A mn . For small n « n c , up to the leading order, we find that the Xq term of Cn(A) is 
exactly n-times of that of Ci(A). From Eq.(15) we easily derive 

Ximp — ~2Sxbuik + subleading order terms, (17) 



where Xbuik ~ T~ 2 m _1 (l/T) is the per-site susceptibility of the bulkEJEil. Very interestingly, the impurity contributes 
a negative susceptibility, which indicates a novel Kondo diamagnetic effect. That means the Kondo coupling dominates 
always over the "molecular field" generated by the bulk ferromagnetic fluctuations. Notice that Eq.(17) is only the 
contribution of the bare impurity. If we take the screening cloud (25 bulk spins which form the bound state with the 
impurity) into account, we find that the total susceptibility of the local singlet is exactly canceled in the leading order. 
That means the polarization effect of the local bound state only occurs in some subleading order, which indicates a 
strong coupling fixed point J* = oo. In fact, the local singlet is much more insensitive to a small external magnetic 
field as we discussed for the ground state. When T — > 0, its susceptibility must tend to zero due to the bound energy 
as shown in Eq.(8). We note the present method is not reliable to derive the total susceptibility of the local singlet 
but the above picture must be true. The same conclusion can be achieved for arbitrary--/ > 0. 

Specific heat of the local composite. In the framework of the local Fermi-liquid theoryEj, the Kondo effect is nothing 
but the scattering effect of the rest bulk particles (N — 25) off the local-spin-singlet composite or equivalently, the 
polarization effect of the local composite due to the scattering. Taking the boundary string into account, the BAE of 
the bulk modes can be rewritten as 

(. \ 2(iV — 2S) AI — 2S 

h Z n -MM TT A 3 -Xt-iXj + Xi-i 

Xj + IJ I) X j -X l +iX j + X l +V [i *> 

i<l>M _ X-i(c + S-l)X + i(c-S-l) / A + § \ 4S 

X + i(c + S-l)X-i(c-S-l) \\-%) ' [ ' 

where </>(A) represents the phase shift of a spin wave scattering off the local composite (boundary bound state). When 
5 = 1/2, c — » 1 + + or J — > +oo, 4>(X) — 0. That means one of the bulk spin is completely frozen by the impurity 
and the system is reduced to an N — 1-site ferromagnetic chain. When 5=1/2, l<c<3/2, only Ci(A) is relevant 
and the free energy of the local composite reads 

Fioc = -T I g(X){( 1 (X)-^( 1 (X-ic + t)-^( 1 (X + ic-i)}dX. (20) 



When xq — 0, we have 



Ci(A) - i&(A-tc + i) - ~Ci(A + ic-i) 



3 1 " C 1 

16(c-l) 2 T^- V lnfH rlm"3 ln^(l + m) H . (21) 

v ; 7T ^ L m(m + 2) J v ; v 1 



The sum in the above equation is convergent for large n c . Therefore we can extend it to infinity, which gives the 
low-temperature specific heat of the local composite as 



4 



(22) 



Similar conclusion can be arrived for arbitrary 5 and J > 0. As long as the Kondo coupling is antifcrromagnetic 
(c > 5+1/2), the low-temperature specific heat of the local composite is described by Eq.(22). There is a slightly 
difference between 5=1/2 case and 5 > 1/2 case. For the former when J — * oo, the local singlet is completely frozen 
and Cioc — > 0, while for the later even when J —> oo, Ci oc takes a finite value. This can be understood in a simple 
picture. For 5 > 1, more than one bulk spin will be trapped by the impurity. Even for J — > oo, only one bulk spin (on 
the nearest neighbor site) can be completely frozen and the rest is still polarizable via the bulk fluctuation. We note 
the specific heat of the local singlet is much weaker than that of the Kondo impurity in a conventional metal. This 
still reveals the insensitivity of the local bound state to the thermal Activation. Though the anomalous power law 
Eq.(22) looks very like that obtained in the Luttinger Kondo systemsEm, they are induced by different mechanisms. 
In the present case, this anomaly is mainly due to the strong quantum fluctuation while in the Lutt inger liquid, the 
anomaly is in fact induced by the tunneling effect of the conduction electrons through the impurityBcSnJ. 

For the ferromagnetic coupling case (J < 0), no boundary bound state exists. Even in the ground state, the 
impurity spin is completely polarized by the bulk Sfftos At finite temperature, the critical behavior is locked into 
that of the bulk {C imp ~ Ti, X im P ~ -(T 2 lnT)" 1 ^. 

Similar phenomena may exist in higher dimensions. The antifcrromagnetic Kondo coupling indicates a local poten- 
tial well for the magnons. Therefore, some bound states of the magnons may exist in the ground-state-configuration, 
which indicate the formation of the local spin-singlet. In this sense, the impurity behaves as a diamagnetic center. 
When J < 0, the Kondo coupling provides a repulsive potential to the magnons and no local bound state can exist 
at low energy scales. The impurity must be locked into the bulk. 

In conclusion, we solve the model of a ferromagnetic Heisenberg chain coupled with a boundary impurity with 
arbitrary spin. It is found that as long as the Kondo coupling is antiferromagnetic, (i) the impurity spin behaves 
as a diamagnetic center and is completely screened by 25 bulk spins in the ground state, no matter how large the 
impurity spin is; (ii) The specific heat of the local composite (impurity plus 25 bulk spins which form bound state 
with it) shows a simple power law C; oc ~ Ti. We note that for a finite density of impurities, the local bound states 
are asymptotically extended to an impurity-band of the magnons, which is very similar to that of a ferrimagnetic 
system. The critical behavior may be different from that of the single impurity case. When the impurity density 
rii ~ 1/(25), we expect a spin singlet ground state. 
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